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TECHNICAL  SUMMARY 


Continued  progress  has  been  made  in  the  past  year  in  understanding 
the  unsteady  flow  about  an  airfoil.  Experimental  data  have  been  obtained 
describing  the  unsteady  pressures  on  the  airfoil  surface  for  conditions 
ranging  from  attached  flow  at  zero  angle  of  attack  to  separated  flow  at 
15®  angle  of  attack.  To  complement  these  data,  the  unsteady  flow  field 
that  produces  these  pressures  has  been  measured  for  the  same  range  of 
test  conditions,  but  with  no  airfoil  present  in  the  wind  tunnel.  Boundary 
layer  measurements,  including  ensemble  averages  of  2  velocity  components 
and  3  Reynolds  stresses  have  been  made  over  a  range  of  attached  flow- 
conditions.  Finally,  wake  surveys  have  been  made  for  50  combinations  of 
the  parameters  Reynolds  number,  perturbation  frequency,  angle  of  attack, 
and  downstream  position. 

Interpretation  of  this  considerable  quantity  of  data  is  continuing, 
with  some  progress  having  been  made  in  identifying  the  key  features  of 
the  flow  field,  and  in  determining  the  relative  significance  of  the  numerous 
parameters  involved.  Work  has  continued  on  the  development  of  numerical 
techniques  to  predict  the  unsteady  boundary  layers  and  wakes  based  on  the 
airfoil  pressure  distributions  and  to  predict  the  pressures  based  upon  the 
unsteady  velocity  field. 

The  present  report  is  divided  into  several  largely  independent  sections. 
First,  a  brief  summary  of  the  experimental  apparatus  and  technique  is  given. 
Second,  some  of  the  results  of  the  measurements  of  the  unsteady  velocity  field 
are  presented.  Third,  the  results  of  measurements  made  to  date  on  the 


pressure  field  of  the  stalled  airfoil  are  given.  Fourth,  earlier  work  on 
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the  airfoil  surface  pressures  and  their  relation  to  the  velocity  field  is 
discussed  by  a  copy  of  a  recently  published  paper  (Appendix  A) .  Finally, 
a  copy  of  a  paper  giving  the  current  results  on  the  unsteady  boundary  layer 
is  included  (Appendix  B) .  Since  the  measurements  of  the  wake  were  completed 
only  at  the  end  of  the  year,  a  report  on  this  data  will  be  given  at  a  later 
date. 


EXPER1MENTA.L  APPARATUS  AND  PROCEDURE 


Since  this  experiment  has  been  described  at  some  length  previously 
(1,2,3),  only  a  brief  summary  is  included  here.  As  shown  in  Figure  1, 
the  experimental  apparatus  consists  of  2-dimensional  NACA  0012  airfoil  which 
may  be  mounted  at  fixed  angle  of  attack  between  sidewalls  in  a  low  speed 
wind  tunnel.  Unsteady  flow  is  generated  by  rotating  an  elliptic  cylinder 
located  behind  and  beneath  the  airfoil  trailing  edge.  Test  Reynolds  num¬ 
bers  ranged  from  7  x  10^  to  1.4  x  10^  based  on  airfoil  chord,  w-hile  reduced 
frequencies,  k  =  aic/2Ua.  varied  between  0.5  and  6.4. 

Airfoil  pressures  were  measured  at  35  locations  on  the  upper  and  lower 
surface,  while  either  single  or  cross  hot  wires  measured  local  velocities. 

A  new  computer-controlled  stepping  motor  probe  drive  system  was  introduced 
this  year  for  use  in  taking  wake  profiles.  This  system  both  speeds  data 
acquisition  and  increases  the  maximum  traverse  from  10cm  to  21cm. 

Data  for  both  pressures  and  velocities  are  required  digitally  as 
shown  in  Figure  2  and  described  in  Reference  3.  Data  storage,  processing 
and  display  techniques  are  also  discussed  in  that  reference. 


TESTS  OF  THE  INDUCED  VELOCITY  FIELD 


To  define  the  unsteady  flow  field  generated  by  the  rotating  elliptic 
cylinder,  the  airfoil  is  removed  from  the  wind  tunnel  and  hotwire  anemometer 
measurements  are  made  of  the  vertical  and  horizontal  velocity  components  along 
the  position  of  the  airfoil  chord.  As  reported  in  References  1,  2  and  4 
such  data  were  taken  along  a  line  corresponding  to  an  airfoil  at  0  angle  of 
attack,  at  position  A,  as  shown  in  Figure  3.  Thin  airfoil  theory  was  then 
used  on  the  data  to  predict  the  unsteady  difference  pressures  on  an  airfoil 
at  that  location.  The  prediction  was  then  compared  to  the  measured  pressure 
distribution  on  the  NACA  0012  airfoil  (1,  2,  3,  4.) 

Following  these  initial  measurements,  airfoil  pressure  data  have  also 
been  taken  at  3  additional  positions  (Figure  3).  In  position  B  the  airfoil 
is  rotated  10*  nose  up  about  the  trailing  edge  from  position  A.  In  position 
C  the  distance  between  the  elliptic  cylinder  and  the  airfoil  trailing  edge 
is  reduced  while  maintaining  10°  angle  of  attack.  In  position  D  the  airfoil 
is  rotated  about  the  trailing  edge  to  15*  angle  of  attack,  while  maintaining 
airfoil-ellipse  distance  as  in  position  C. 

Since  the  induced  velocity  field  at  these  additional  positions  differs 
from  that  at  position  A,  an  additional  series  of  elliptic-cylinder  alone 
velocity  measurements  were  required  to  determine  this  difference.  These 
tests  were  performed  in  November  1981,  using  the  digital  data  acquisition 
and  processing  system  that  had  been  developed  since  the  original  analog 
measurements  were  made  (Figure  2,  Reference  3).  The  u  and  v  velocity  com¬ 
ponents  were  measured  at  21  locations  which  covered  the  4  chord  positions. 

Data  were  taken  at  2  Reynolds  numbers  (7  and  10  x  10^  based  on  airfoil  chord) 


and  at  steady  ai  3  5  unsteady  elliptic  cylinder  frequencies  (0.5<k<6.4). 


Allowing  for  repeated  runs,  26S  test  points  in  (position,  Re,  k)  space 
were  taken. 

The  qualitative  features  of  the  ensetible  average  induced  upwash  (v) 
velocity  field  may  be  seen  in  Figures  4  and  5.  These  figures  are  3-dimen¬ 
sional  plots  of  the  upwash  as  a  functio.n  of  time  and  distance  along  the 
position  of  the  airfoil  chord.  Time  has  been  nondimensionalized  so  that 
720°  corresponds  to  1  elliptic  cylinder  revolution.  Distance  from  the 
leading  edge  position  is  nondimensionalized  by  the  airfoil  chord,  and 
velocity  by  the  freestream  velocity.  Both  figures  are  at  Re  =  7  x  1C'' 
and  chord  position  A.  Figure  4  is  for  k  =  1.0,  and  shows  the  typical 
character  of  the  upwash:  both  mean  and  unsteady  velocities  increase 
sharply  as  the  trailing  edge  position  and  the  elliptic  cylinder  are 
approaclied;  the  waveform  is  roughly  sinusoidal  in  time;  and  the  shift 

in  phase  with  distance  along  the  chord  is  small.  Figure  5,  for  k  =  3.9, 
is  qualitatively  similar,  but  is  characterized  by  higher  mean  upwash  velocities 
and  lower  unsteady  amplitudes. 

Mean  upwash  distributions  for  the  two  extreme  steady  elliptic  cylinder 
orientations  at  position  C  are  shown  in  Figure  6.  The  upwash  for  the  horizontal 
ellipse  shown  is  based  on  a  potential  flow  calculation.  This  upwash  is  used 
as  the  reference  upwash,  and  is  needed  for  the  cross  hot-wire  anemometer 
calibration  (1)  .  The  upwash  for  the  vertical  ellipse  shown  is  a  sum  of  t.his 
reference  upwash  and  the  measured  difference  in  upwash  between  the  vertical 
and  horizontal  cases.  The  dashed  line  in  Figure  6  is  a  potential  flow  calcu¬ 
lation  for  the  vertical  ellipse.  Naturally,  it  does  not  account  for  the 
large  regions  of  separated  flow’  which  exist  behind  the  vertical  ellipse. 

The  procedure  of  adding  a  calculated  reference  upwash  to  a  measured  differ¬ 
ence  in  upwash  was  followed  for  all  steady  and  unsteady  cases  and  is  explained 
in  Ref.  1.  Note  that  this  procedure  includes  tunnel  non-uniformities  in 


the  reference  flow. 


To  get  quantitative  information  on  the  upwash  distributions,  the 
ensemble  average  waveform  such  as  in  Figure  4  and  5  is  Fourier  analyzed, 
and  the  amplitude  and  phase  lag  of  the  harmonics  of  the  ellipse  rotation 
rate  are  found.  The  fundamental  harmonic  is  defined  to  be  twice  the 

C 

ellipse  rotation  rate.  Figure  7  presents  amplitude  data  for  Re  =  7  x  1C  , 
k  =  1,  position  C.  Haraionic  0  is  the  mean,  0.5  the  subharmonic,  and  1 
the  fu.ndamental .  The  distributions  of  the  mean  and  of  the  fundamental 
confirm  the  qualitative  features  discussed  above,  while  those  for  the 
sub  and  higher  harmonics  confirm,  the  essentially  sinusoidal  nature  of  the 
flow.  Note  that  while  the  original  prediction  of  unsteady  airfoil  pressure 
(Ref.  1)  considered  only  the  first  harmonic,  the  effect  of  the  other 
harmonics  may  be  added  if  desired. 

Figure  8  presents  the  same  quantities  for  k  =  3.9.  The  major  differ¬ 
ences  from  the  k  =  1  results  in  Figure  7  are  the  reduced  unsteady  amplitudes, 
the  somewhat  larger  mean  amplitudes,  and  the  smaller  mean  amplitude  for 
x/c  >  1. 

Phase  lag  distributions  for  position  C,  Re  =  7  x  10^  and  0.5  <  k  <  6.4 
are  shown  in  Figure  9  as  a  function  of  x/c.  A  phase  lag  of  0  means  that 
the  minimum  upwash  occurs  when  the  elliptic  cylinder  is  horizontal.  All 
of  the  frequencies  are  characterized  by  nearly  constant  phase  lag  for  x/c  <  0.8, 
with  an  increase  in  phase  lag  for  x/c  >  0.8.  For  k  <  1,  this  increase  is 
approximately  20°,  where,  as  before,  720°  corresponds  to  one  elliptic  cylin¬ 
der  rotation.  At  k  =  2  the  increase  is  of  the  same  magnitude,  but  the  phase 
lag  is  about  40°  larger  at  all  points  along  the  chord  than  at  k  <  2. 

For  k  >  2  the  phase  lag  at  all  points  drops  again,  with  values 


for  k  =  6.4  beinc  appro.xirnately  60'  loss  tha.r,  th'--  values  for  k  =  1  at  the 
same  position.  The  increase  in  ph.ase  between  x/c  =  0.6  a.nd  >;,'c  =  1.1  is 
larger  when  k  >  2,  reaching  45'  for  .k  =  6.4. 

For  the  present  purposes,  the  details  of  the  flow  field  are  less 
important  than  the  fact  that  the  phase  is  relatively  constant  along  a 
chord  line,  for  a  given  frequency,  Rey.nolds  number  and  trailing  edge  posi¬ 
tion.  In  contrast,  if  the  flow  field  was  due  to  a  sinusoidal  gust  convected 
at  a  uniform  velocity  U  ,  the  phase  distribution  would  be  described  by: 


^  ^  d: _ 2k 

dx  U  '  d  (x/c)  (U  /u^) 

9  9 

which  yields  a  phase  difference  between  the  leading  and  trailing  edges 

of  lid  =  2kU  /U  radians.  For  U  =13  and  k  =  .5,  di  =  57.3',  while 
oc  c  oc 

for  k  =  6.4,  Ai  =  734®,  almost  fifteen  times  larger  tha.n  the  measured 
values.  Therefore  the  assumption  of  a  fixed  or  at  best  a  slowly  converted 
disturbance  see.mis  reasonable. 

Note  that  a  correction  to  the  unsteady  thin  aiifoil  theory  to 
account  for  the  small  convective  component  of  the  upwash  may  be  made  (4) . 

This  correction  makes  an  insignificant  change  in  the  resulting  pressure 
distributions  under  the  conditions  studied  here. 

To  obtain  a  more  global  picture  of  the  dependence  of  the  induced 
upwash  on  the  parameters  Re,  a,  k,  and  position,  the  upwash  at  the  trailing 
edge  location  was  selected  as  a  reference.  This  choice  reflects  both  the 
importance  of  the  trailing  edge  in  determining  airfoil  flows,  and  the 
fact  that  the  maximum  upwash  occurs  there.  Figure  10  shows  the  mean 
flow  angle  (arctangent  (v/u))  with  respect  to  the  airfoil  chord  line,  plotted 


against  reduced  frequency  for  the  4  positions  and  the  2  Reynolds  numbers  studied. 


In  general  the  mean  flow  angle  at  the  trailing  edge  shows  little  cerjendcncc- 
on  Re  or  k,  being  influenced  chiefly  by  the  airfoil  angle  of  attack  and  the 
distance  to  the  elliptic  cylinder.  Thus  results  for  position  B  are  lO"^  higher 
than  for  A  due  to  angle  of  attack,  and  results  for  C  are  5-7°  higher  than  for 
D  due  to  reduced  trailing  edge-elliptic  cylinder  distance. 

Figure  11  shows  the  amplitude  of  the  fundamental  harmo.nic  of  the 
flow  angle  for  the  same  conditions.  The  Reynolds  .number  dependence  is 
larger  than  that  seen  in  the  mean  angle,  but  is  still  small  compared  to 
the  frequency  depende.nce.  The  Reynolds  number  effect  is  larger  nearer 
to  the  elliptic  cylinder,  where  the  details  of  the  flow  over  it  are  more 
important  than  they  are  further  away.  Due  to  the  increased  acceleration 
needed  to  move  the  fluid  at  higher  frequencies,  the  amplitude  is  reduced. 

The  dependence  on  position  is  strong,  witn  amplitudes  for  positions  C 
and  D  50%  larger  than  for  A  and  B  at  k  <  1,  and  up  to  200%  larger  at 
k  =  6.4.  The  p.hase  lag  at  the  trailing  edge,  as  shown  in  Figure  12, 
has  a  similar  weak  dependence  or.  Reynolds  number,  and  a  strong  dependence 
on  frequency  and  distance.  Phase  lags  appear  to  approach  0  as  k  approaches 
0,  in  agreement  with  steady  incompressible  flow,  where  the  propagation 
speed  is  infinite. 

To  conclude,  these  more  complete  and  detailed  studies  of  the  induced 
velocity  fields  largely  confirm  trends  observ'ed  earlier,  extending  the.m  to  the 
wider  range  of  parcimeters  (k,  a,  Re,  position)  currently  studied.  The 
field  depends  smoothly  on  space  and  time,  with  maximum  amplitudes 
concentrated  near  the  position  of  the  airfoil  trailing  edge,  and  unsteady 
components  dominated  by  the  fundamental  harmonic.  Reynolds  number  dependence 


is  small  away  from  the  immediate  vicinity  of  the  elliptic  cylinder. 


The  freq-ceriCy  dependence  of  the  unsteady  corriponent  of  the  flow  is  always 
significant.  No  qualitative  changes  were  seen  in  the  field  due  to  changes 
in  airfoil  position  relative  to  the  ellipse.  Cnsteady  a-T^plinudes  were, 
as  desired,  significantly  increased  by  moving  the  trailing  edge  closer  to 
the  elliptic  cylinder.  Phase  lag  distribution  along  the  chord  was  studied 
in  more  detail  than  previously  possible.  The  chordwise  variation  of  phase 
was  found  to  be  small  and  is  much  smaller  than  that  expected  if  a  distur¬ 
bance  is  convected  at  the  free  stream  velocity. 


PRESSURE  DISTRI&'JTICr;  ON  THE  STALLED  MRTOIL 

In  uniform  steady  flow,  the  NACA  0012  airfoil  has  been  found  to  stall 
at  an  angle  of  attack  of  approximately  16°,  at  a  lift  coefficient  of  1.5  (7). 
These  values  are  for  Re  =  3  x  10  ,  the  lowest  value  given.  In  the  present 
study  the  airfoil  is  in  the  non-uniform  unsteady  flow  generated  by  the 
rotating  elliptic  cylinder,  so  that  the  conditions  causing  stall  are 
different,  as  will  be  discussed  below. 

Stall  was  observed  for  position  C  at  a  =  10°  only  for  k  =  6.4  and 
Re  =  7  X  10^,  and  for  all  frequencies  and  Reynolds  n^u-tbers  studied  at 
position  D  at  a  =  15°.  The  stall  was  characterized  by  a  separation  of  the 
boundary  layer  on  the  upper  surface  starting  near  the  leading  edge  and  ex¬ 
tending  over  the  entire  surface.  No  dyna.mic  stall,  or  oscillatory  separation 
and  reattachment  was  seen.  The  unsteady  pressures  generated  by  the  rotating 
elliptic  cylinder  propagated  through  the  separated  flow  region,  but  the 
essential  character  of  the  flow  remained  constant. 

The  first  case  to  be  discussed  is  the  separated  flow  that  may  occur 
at  a  =  10°,  k  =  6.4  and  Re  =  7  x  10^.  Under  these  conditions  both  attached 
and  separated  flow  states  are  possible.  The  mean  pressure  distributions  for 
these  two  states  are  shown  in  Figure  13.  The  attached  flow  is  characterized 

by  a  suction  peak  on  the  upper  surface  at  x/c  =  .025  of  c  =  -6.2,  while  the 

P 

separated  flow  has  a  nearly  uniform  mean  upper  surface  pressure  of  -1  <  c  <  -0. 

P 

The  lower  surface  pressures  are  almost  the  same  for  the  two  states,  largely 
unaffected  by  the  stall.  The  flow  may  remain  attached  for  times  up  to  an  hour, 
or  until  a  disturbance  large  enough  to  trip  the  flow  passes  the  airfoil,  Once 
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secsrated  the  flow  has  never  been  observed  to  spontanecusl y  reattacr. •  If 
the  reduced  frequency  is  reduced  to  approximately  k  =  the  flow  is  re¬ 
attached  and  the  frequency  may  aqain  be  increased  tc  k  =  6.4.  Tnis 
hysteresis  seems  typical  of  stalled  flows.  (Ref.  8). 

Ensemble  averaged  unsteady  pressure  distributions  may  be  compared  to 
the  2  cases.  Approximately  400  cycles  of  averaging  are  required  tc  define 
the  flow  in  the  separated  case,  while  only  200  are  sufficient  for  the 
attached  flow.  Figure  14  shows  upper  and  lower  surface  fundamental  harmonic 
amplitude  distributions.  The  lower  surface  distribution  is  seen  to  be 
relatively  unaffected  by  the  separation,  with  pressure  amplitudes  u.niformly 
increased  by  about  15%  by  separation.  The  upper  surface  distribution  is 
quite  different,  as  separation  appears  to  superimpose  a  damped  sinusoidal 
oscillation  on  the  smoothly  varying  attached  distribution.  The  phase  lag 
distributions  in  Figure  15  show  similar  behavior.  The  lower  surface  is 
virtually  unaffected,  while  an  oscillation  in  phase  is  superimposed  on  the 
upper  surface  for  x/c  <  0.8, 

By  inserting  phase  jumps  of  2f  at  x/c  =  0.05  and  0.3  a  region  of  monotonic 
increase  in  phase  lag  for  .025  <  x/c  <  0.4  is  obtained,  such  as  that  shown 
later  in  Figure  24.  A  linear  fit  in  this  region  gives  results  similar  to  those 
of  a  disturbance  propagating  downstream  at  a  phase  velocity  of  0.5  U^.  That 
the  situation  is  more  complex  than  a  simple  propagating  wave  is  apparent  from 
the  amplitude  distribution,  but  the  picture  of  pressure  disturbances  from  the 
lower  surface  moving  around  the  leading  edge  and  propagating  downstream  after 
the  stagnation  point  may  be  useful.  A  similar  result  was  obtained  by  a  study 
of  the  fluctuating  pressures  on  a  steady  airfoil  (8). 

At  airfoil  position  D,  where  the  angle  of  attack  is  15°,  all  5  unsteady 
frequencies  (0,5  <  k  <  6.4)  studied,  plus  the  case  of  a  steady  vertical 


elliptic  cylinder  at  both  test  Reynolds  numbers  of  7  and  ID  x  10“  were  stalled. 
Only  the  case  of  a  steady  horizontal  elliptic  cylinder  had  attached  flow. 

Figure  16  shows  the  mean  pressure  distributions  at  Re  =  10^  for  all  studied 
frequencies.  As  in  the  previous  stalled  case,  the  lower  surface  is  relatively 
unaffected,  while  the  upper  surface  pressure  distribution  is  flat  at  -l<c^<  -0.5. 
Frequency  dependence  is  exhibited  pri.marily  on  the  upper  surface:  leading 
edge  pressures  are  higher  for  k  <  2,  a  low  pressure  bump  develops  at  k  =  2, 

and  more  variation  with  x  is  seen  near  the  leading  edge  for  k  >  2. 

A  more  quantitative  picture  is  obtained  by  examining  the  mean  and 
unsteady  lift  coefficients.  Figure  17  presents  the  mean  lift  coefficient  as 
a  function  of  frequency  for  both  the  stalled  a  =  15®  data  and  the  attached 
a  =  10°  data  for  the  same  trailing  edge  position  relative  to  the  ellipse. 

The  attached  flow  mean  lift  coefficie.nt  increases  with  k,  primarily  due 
to  the  increased  circulation  produced  by  spinning  the  elliptic  cylinder 
faster,  reading  a  maximum  of  C  =  1.46  for  k  >  5.  This  value  is 

more  or  less  equal  to  the  maximum  C.  for  the  NACA  0012  airfoil  in  steady 

Li 

uniform  flow  (7).  However,  if  the  unsteady  lift  coefficient  is  added  to  the 

mean,  the  instantaneous  values  of  C  become  as  large  as  1.6,  greater  than  the 

steady  maximum.  These  results  may  be  compared  with  those  obtained  by  rotating 

the  NACA  0012  in  uniform  flow,  such  as  by  McCroskey,  Carr  and  McAlister  (9). 

6 

In  these  tests  the  Reynolds  number  is  somewhat  larger.  Re  =  2.5  x  10  ,  the 
freouency  lower,  k  ^  0,25,  and  the  amplitude  much  greater,  5  ^  ^  ^  ■ 

Instantaneous  maximum  lift  coefficients  were  found  to  be  much  greater  than  the 
maximum  in  steady  flow,  but  when  averaged  over  the  airfoil  rotation  period, 
the  mean  lift  coefficient  was  much  closer  to  the  steady  flow  value. 

Mean  C  increases  with  Reynolds  number,  also  in  agreement  with  the  high  lift 
results  in  Ref.  7.  The  stalled  flow  lift  coefficients  are  always  lower,  with  the 
difference  between  being  less  than  10%  for  k  ^  2,  and  becoming  30%  lower  for  the 
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deeper  stall  case  at  k  >  5.  As  in  the  pressure  distributions  in  Figure  16, 
the  behavior  changes  significantly  at  k  =  2.  The  reasons  for  this  change 
are  not  known. 

The  unsteady  lift  coefficient  is  relatively  unaffected  by  the  stall  of 
the  mean  flow.  Figure  18  shows  the  amplitude  of  the  fundamental  harmonic 
for  the  cases  discussed  above.  The  attached  flow  data  re  in  qualitative 
agreement  with  earlier  results  for  positions  A  and  B  (Ref.  3),  and  after 
the  mean  flow  separates  the  first  harmonic  still  follows  the  same  trends: 

C  drops  for  k  >  1,  has  a  minimum  at  k  ^  2  and  then  increases  for  k  >  2. 

Li 

Phase  lags  for  C  ,  as  shown  in  Figure  19,  are  also  qualitatively  similar 

Ij 

for  the  attached  and  separated  mean  flows,  with  the  prLmary  difference 
being  a  45°  higher  phase  lag  for  k  <  1. 

The  ensemble  averaged  pressure  distribution  is  shown  in  the  next  series 
of  figures,  which  are  3-dimensional  plots  of  the  pressure  coefficient  as 
a  function  of  non-dimensional  time  and  distance  from  the  airfoil  leading 
edge.  Each  curve  is  an  average  of  400  cycles.  Figure  20  shows  the  pressures 
on  the  lower  surface  for  k  =  1,  a  =  15°,  and  Re  =  10^.  As  would  be  expected 
from  the  curves  discussed  earlier,  the  pressure  distribution  is  similar  to 
lower  surface  distributions  for  unstalled  cases,  and  is  characterized  by  large 
unsteady  amplitudes  and  small  phase  differences  along  the  chord  (Ref.  3).  Figure 
21  shows  the  upper  surface  for  the  same  conditions.  Note  that  the  x/c  axis 
has  been  reversed  for  clarity  in  Figure  21.  In  the  trailing  edge  region  the  pres¬ 
sure  is  similar  to  that  for  unstalled  flows,  with  large  amplitudes  and  small  phase 
shifts.  From  the  leading  edge  to  midchord  the  amplitude  varies  relatively  rapidly, 
and  the  phase  lag  increases  monotonically .  This  is  consistent  with  the  results 
shovm  in  Figure  14  and  15  at  a  =  10°.  The  difference  pressure  coefficient, 

C  ,  C  ,  is  shown  in  Figure  22.  It  includes  both  the  nearly  constant 

phase  lower  surface  wave  form  and  the  rapidly  varying  upper  surface  disturbance. 


The 


resulting  in  a  pattern  quite  different  froni  that  of  an  unstalled  flow, 
interface  between  these  two  regions  near  ir.idchord  appears  for  this  reason 
to  have  a  frequency  doubling. 

Data  for  higher  reduced  frequencies  are  sLn-.ilar,  but  not  as  clearly  defined, 
due  primarily  to  the  more  rapid  phase  variation,  a  result  of  a  decrease  wave¬ 
length  to  chord  ratio.  No  significant  Reynolds  number  effect  was  seen  on 
the  ensemble  averaged  pressure  coefficients. 

The  exception  to  this  behavior  is  the  region  near  the  trailing  edge 
of  the  upper  surface  at  k  =  0.5,  as  seen  in  Figure  23  for  Re  =  7  x  10^. 

A  frequency  doubling  develops  for  x/c  >  0.9,  due  perhaps  to  a  coupling  with 

the  vortex  shedding  from  the  slowly  rotating  elliptic  cylinder.  The  same 

6 

effect  IS  present  at  Re  =  10  . 

As  discussed  above  for  the  cx  =  10®  data,  one  interpretation  of  the 
pressure  distribution  on  the  upper  surface  is  in  terms  of  the  propagation 
of  a  travelling  disturbance  downstream  from  the  leading  edge.  The  ensemble 
average  pressure  may  be  represented  as 

-  n,  ,0. 

<p>  (x,t)  =  p{x)  +  p^  (x ) cos  (cot  -  d(x))  +  P2(x,t), 


where  contains  harmonics  other  than  the  fundamental.  For  a  disturbance 

propagated  at  a  constant  speed  U  ,  the  phase  lag  4'(x)  is  cox/U  .  There- 
^  ^  ^  ^  wave  wave 

fore  the  wave  velocity  is  defined  as  U  =  a)/(d6/dx),  or  in  non-dimensional 

•'  wave 

terms  U  /U_  =  2k/ (cd4>/dx) .  If  the  mean  flow  velocity  U  (x)  is  introduced, 
wave  00  '  g 

its  ratio  with  this  wave  velocity  becomes 


U  /U 
wave  e 


2k 
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Figure  24  shows  the  phase  lag,  C(x),  for  two  representative  cases  at 
0 

Ct  =  15°,  Re  =  10  .  For  k  =  1  the  phase  increases  smoothly  and  moriotcnically , 
but  not  linearly  as  in  the  simple  model  above.  For  k  =  4  the  data  is  similar, 
but  more  complex,  as  seen  in  the  jump  in  phase  near  x/c  =  .5  and  in  the  ma.xi- 
mum  at  x/c  =  .75,  followed  by  a  slight  drop  in  phase  lag  as  the  trailing 
edge  is  approached. 

Figure  25  shows  the  linear  propaqation  velocity  U  /U  for  these  same 

wave  e 

two  cases.  For  k  =  1  the  variations  are  smooth,  U  is  always  positive 

wave 

and  normally  in  the  range  0.2<U  /U  <1.2,  consistent  with  a  disturbance 

wave  e 

propagating  downstream.  The  situation  is  much  more  complex  at  k  =  4,  where 
the  wave  speed  varies  much  more  rapidly,  reaching  extremes  of  +6U^  and 
The  negative  wave  velocities  for  x/c  >  0.85  seerr.  to  indicate  an  upstream 
propagation  in  this  region.  The  interaction  between  this  disturbance  and 
that  propagating  downstream  from  the  leading  edge  may  account  for  the  unusual 
behavior  for  0.6  <  x/c  <  0.8.  Another  possible  cause  is  interaction  with 
the  two-dimensional  channel  surrounding  the  airfoil. 

The  propagation  of  disturbances  in  separated  flow  is  obviously  a  very 
complex  process,  and  not  a  simple  travelling  wave  or  sum  of  travelling  waves, 
as  with  a  convected  gust.  The  process  is  highly  frequency  dependent,  with 
the  complexity  increasing  rapidly  for  k  >  1. 

In  conclusion,  some  of  the  main  characteristics  of  the  separated  flow 
about  an  airfoil  in  our  particular  unsteady,  non-uniform  flow  field  have 
been  discussed.  Many  of  these  characteristics,  in  particular  those  dealing 
with  conditions  for  separation,  the  mean  pressure  distribution,  and 
the  behavior  of  the  lift  coefficient  appear  to  agree  with  results  of  earlier 
attached  flow  data  and  other  results  for  stalled  flow.  The  unsteady  pressure 
distribution  at  the  fundamental  harmonic  seems  relatively  unaffected  by 
the  outer  mean  flow  separation.  However,  other  characteristics,  primarily 
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those  concerned  with  the  propagation  of  unsteady  pressures  in  the  separated 
flow  are  more  complex  and  resist  straightforward  explanations.  More 
effort,  both  conceptual,  in  interpreting  the  existing  data,  and  experimental , 
primarily  in  adding  velocity  measurements  to  go  with  the  pressure  data, 
is  necessary  in  order  to  reach  an  improved  understanding  of  the  process (es). 
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LIST  OF  SYMBOLS 

airfoil  chord 

2 

section  lift  coefficient,  L/(l/2pU^C) 

2 

pressure  coefficient,  (P-P  ^  )/(1/2dU^) 

^  static  ^  “ 

reduced  frequency,  tic/2U^ 

section  lift 

pressure 

freestream  static  pressure 
time  averaged  pressure 
ensemble  averaged  pressure 
fundamental  harmonic  of  <P> 
other  harmonics  of  <F> 

Reynolds  number , 
time 

tangential  velocity 

external  velocity,  at  boundary  layer  edge 

gust  velocity,  propagation  velocity  of  a  disturbance 

normal  velocity 

tangential  coordinate 

normal  coordinate 

angle  of  attack,  with  respect  to  freestream  velocity 
kinematic  viscosity 
density 
phase  angle 


radian  frequency 


Status  of  Research 


1.  Earlier  conclusions  as  to  the  nature  and  superpcsability  or  the 
unsteady  pressure  field  on  the  mean  pressure  field  and  the  applicability 
of  unsteady  thin  airfoil  theory  have  been  confirmed  and  extended  from 
the  zero  angle  of  attack  to  higher  angles  of  attack. 

2.  A  series  of  tests  on  the  pressure  distribution  in  stalled  flow  have 
resulted  in  a  collection  of  data  revealing  many  characteristics 

of  the  flow.  Further  thought,  interpretation  and  investigation  is 
required,  in  particular  with  regard  to  the  origin  and  nature  of  pressure 
disturbances  in  this  flow. 

3.  Boundary  layer  measurements  have  been  taken  over  a  wide  range  of 
pressure  gradients  and  frequencies,  allowing  determination  of  some 

of  the  key  features  of  the  unsteady  boundary  layer  and  of  the  influence 
of  the  important  parameters.  More  work  is  req'uired  in  order  to  measure 
profiles  at  more  stations  at  each  configuration,  to  get  a  better  picture 
of  the  development  of  the  boundary  layer.  Continued  development  of 
numerical  techniques  to  predict  the  boundary  layer  with  reasonable 
accuracy  is  also  required. 

4.  Wake  profiles  have  been  taken  over  a  similarly  wide  segment  of  parameter 
space,  and  will  provide  a  data  base  for  increasing  understanding  of  this 


little  investigated  portion  of  the  unsteady  flow. 


Development  of  a  shear  gauge  to  measure  unsteady-  skin  friction  has 
been  slowed  by  the  persistent  drift  in  the  electronic  circuits. 
Presently  this  effort  is  on  hold  due  to  lack  of  a  suitable  research 
assistant.  We  are  considering  the  use  of  other  techniques  that  have 
been  developed  recently  to  determine  this  important  parameter. 
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The  following  people  have  contributed  to  this  project  in  the 
past  year : 


1.  Professor  E.E.  Covert 

2.  P.F.  Lorber 

3.  R.  Lee  (up  to  June  19S2) 

4.  C.M.  Vaczy  (after  June  1982) 


Principal  Investigator 
Research  Assistant 
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Interactions 


The  paper  "Unsteady  Turbulent  Boundary  Layers  in  Adverse  Pressure 
Gradients"  was  presented  at  the  AIAA/A5.ME  3rd  Joint  Theririophy  sics ,  Fluids, 
Plasma  and  Keat  Transfer  Conference,  held  Ju,ne  7-11,  19S2  in  St.  Louis, 

MO. 

We  have  discussed  the  unsteady  data  with  Professor  Dennis  Whitehead 
of  the  WTiittle  Laboratory  with  a  view  of  application  of  his  Finite  Element 
Technique  to  assist  interpretation  of  data. 


1)  Test  Section 

2)  NACA  0012  Airfoil 

3)  Rotating  Elliptic  Cylinder 

4)  Drive  Motor  (0-3300  rpm) 

5)  2-D  Sidewalls 

6)  Pitot-Static  Probe 
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REDUCED  FREQUENCY 
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